Dipole oscillations of confined lattice bosons in one dimension 
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We study the dynamics of a non-integrable system comprising interacting cold bosons trapped in 
an optical lattice in one-dimension by means of exact time-dependent numerical DMRG techniques. 
Particles are confined by a parabolic potential, and dipole oscillations are induced by displacing the 
trap center of a few lattice sites. Depending on the system parameters this motion can vary from 
undamped to overdamped. We study the dipole oscillations as a function of the lattice displacement, 
the particle density and the strength of interparticle interactions. These results explain the recent 
experiment C. D. Fertig et al, Phys. Rev. Lett. 94, 120403 (2005). 



Recent experiments with cold atoms P, 0, Q have 
provided realizations of non-equilibrium quantum many- 
body systems, allowing to address a number of funda- 
mental questions. For example, the integrability of a 
many-body system has been demonstrated in Ref. 0, 
via the inhibition of thermalization in a one-dimensional 
Bose gas, which opened the way to theoretical studies of 
the relaxation dynamics of non-equilibrium many-body 
systems [HI . The dynamics of non-integrable systems has 
been recently explored experimentally in Refs. [1,01 using 
interacting cold bosonic atoms trapped in an array of one- 
dimensional optical lattices and confined by a parabolic 
potential. Dipole oscillations were induced by displac- 
ing the center of the parabolic potential, and the dipole 
dynamics was studied by monitoring the position of the 
center of mass. A sudden transition from a regime of 
undamped motion to a regime of strongly damped mo- 
tion was observed on increasing the lattice depth. Since 
damping of the center of mass oscillations is due to exci- 
tations in the optical lattice, the results obtained in 0, Q 
have provided precious diagnostic of the dynamical cor- 
relations of the many-body system, and thus have stim- 
ulated considerable theoretical interest 0, S, 0| ■ 

Good agreement with the experimental results in Q 
has been obtained in the regimes of very weak and very 
strong interactions 3, where mean-field and extended 
fermionization techniques apply. However, it remains 
a fundamental challenge to understand the dipole dy- 
namics in the regime of intermediate interactions, where 
the sudden localization transition occurs and the sub- 
tleties of one-dimensional (ID) correlations do not allow 
(semi-) analytical treatments. With the aim to provide a 
coniprehensive explanation of the experiment of Fertig et 
al. [J|, in this letter we study the dipole oscillations by 
means of a numerically exact time-dependent density- 
matrix-renormalization-group technique (tDMRG), see 
also We find very good agreement with the ex- 

perimental results in the interesting regime of inter- 
mediate interactions. These results demonstrate that 



time-dependent numerical simulations with tDMRG have 
reached the same accuracy of current experiments with 
cold gases in the strongly correlated regime and thus rep- 
resent a unique theoretical tool for quantitative compar- 
isons and predictions for experiments in the cold atoms 
context. 

The experiment in was performed in a parame- 
ter regime where the use of the following Bose-Hubbard 
Hamiltonian is microscopically justified 
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The first term on the r.h.s. of Eq.(IT]) describes the tun- 
neling of bosons between neighboring sites with rate J {j 
labels the sites on the lattice). The second term is the 
parabolic potential with curvature 51; S{t) is a sudden dis- 
placement of the trap center, do{t) = S Q{t) (with 8(i) 
the Heaviside function), and nj = b'jbj is the density op- 
erator with bosonic creation (annihilation) operators bj 
(hi). The last term is the onsite contact interaction with 
energy U [llj, (we set h= 1). 

The sudden displacement on the trap center causes 
dipole oscillations of the bosons which can be analyzed 
experimentally by monitoring the time evolution of the 
Center Of Mass (COM) Xcom = J2jj{nj)/N, with N 
the number of particles. The experiment of Ref. 0] was 
performed on a array of one-dimensional optical lattices 
where the number of particles in each ID lattice var- 
ied from ~ 80 to zero. Thus, in order to provide a 
comprehensive and quantitative comparison with the ex- 
perimental data, here we analyze the dipole dynamics as 
a function of 6, U/J, and the number of bosons N. We 
find that overdamped motion can occur as a function of S 
for arbitrarily small interactions. Fig. [21 while in general 
sizeable interactions tend to extend the parameter region 
where localization occurs [12| . For a given Q/ J damping 



2 




FIG. 1: Relevant density distributions [panels (a-c)], see text, 
and snapshots of the corresponding dipole dynamics [panels 
(d-f)]. (a-b) Density distribution for A'' = 11 and 15 particles, 
respectively, for fl/J = 0.05623. In each panel, the dashed 
and solid lines are U/J = 1 and 20, respectively. The solid 
line in panel (b) corresponds to a Mott insulator, (c) Den- 
sity distribution for N ^ 23, Q/J = 0.4, and U/J = 1 and 
20 (dashed and solid lines, respectively). The solid line corre- 
sponds to a cake- like structure, (d-f) Snapshots of the density 
distribution for the cases (a-c), with U/J = 20 and 5 = 4, at 
times tJ = 0, 30, 40 and 50. The dynamics of a few atoms 
in the Mott and cake-like configurations is frozen, however, 
residual oscillations can persist in the latter, see text. 



is found to depend exponentially on U/J, and to be fa- 
vored for small N. Figure [HJa), where the damping rate 
is shown as a function of the interaction and the num- 
ber of bosons and, most important. Fig. [3l^b), where we 
compare our numerics with the experimental data finding 
very good agreement in the intermediate range of inter- 
actions, allow for a new explanation of the experiment of 
Ref. based on the role of lattices with different N. 

Three regimes are of interest for the dipole dynamics 
[sec Fig. [1]: a) for 4 J > n{N/2)'^ the density distribu- 
tion is Gaussian or Thomas-Fermi-likc for 4J ^ U and 
4J ~ U, respectively, and for U ^ onsite densities 
are smaller than one; b) for U > il{N/2)'^ > 4J a Mott 
insulator with one particle per site is formed at the trap 
center; c) for Q{N/2)^ > U > AJ a sheU structure is 
formed with a density 1 < n-,- < 2 at the trap center, sur- 
rounded by a Mott-insulator with one particle per site. 
All the situations above occur in the experiment, since 
N varies from one lattice to another. Therefore, in the 
following wc arc first interested on the dynamics of model 
systems as those in Fig. [1] which exemplify all three cases 
a), b) and c) above while still allowing for an extensive 
analysis in terms of all parameters N, Q/ J and U /J, and 
then we address the experiment of Ref. in the most 
interesting regime U/ J > 4. 

The results presented below have been obtained by 
means of a tDMRG algorithm with a second order Trot- 
ter expansion of H, and time-steps O.OIJ Q. We take 
advantage of the conserved total number of particles N 



projecting on the corresponding subspace; the truncated 
Hilbert space dimension is up to to = 100, while the al- 
lowed number of particles per site is D — 5. All results 
below are found to be independent of this choice. 

We first focus on the dipole dynamics as a function of 
the trap displacement S, in the regime of weak interac- 
tions. In this regime, mean-field theory predicts a sud- 
den transition between undamped and overdamped mo- 
tion via a dynamical instability at a critical displacement 
Sc — \f%jj^ [l3|. This value for 5c can be understood 
by employing the exact solution of Eq. ([T]) in the non- 
interacting limit [l3|. For energies E < AJ the single- 
particle eigenstates of H{t ~ 0) are harmonic-oscillator- 
like modes extended around the center of the parabolic 
trap. However, for E > AJ particles are Bragg-scattered 
by the lattice, and perform Bloch-like-oscillations cen- 
tered far from the trap center IJj. The particle lo- 



calization corresponds to the population of these latter 
high-energy modes, which becomes significant for dis- 
placements 5 >j 5c, [S]. Our numerical results in the 
limit of weak interactions are shown in Fig. [2^ a-c), where 
dipole oscillations of the center of mass Xcom are shown 
as a function of time t, for different values of the dis- 
placement 5. In the simulations, as initial condition we 
use the ground-state wavefunction of the undisplaced po- 
tential, shifted by 5 lattice sites. On increasing 5, the 
dynamics changes from undamped to damped, and the 
particles oscillate around the trap center. On increas- 
ing further the displacement [5 > 5 in panels (a-b)] the 
oscillations are overdamped, and the COM slowly drifts 
towards the trap center or clings to the borders of the 
trap [case with iV = 23 of panel c)]. This behavior corre- 
sponds to the localization transition predicted by mean- 
field theory. However, Fig. [2] shows that quantum fiuc- 
tuations, properly accounted for by the tDMRG, smear 
out the transition into a smooth crossover between the 
undamped and the overdamped regimes. 

Having established a connection with known results 
in the mean-field regime, we now present exact results 
for the particle localization in the interesting case of 
stronger interactions C// J > 1 and 5 < 5c- We first focus 
on model systems and fix (5c = 6 and the displacement 
5 = 1 < 5c, such that for small interactions U/J < I 
the dynamical instability discussed above does not occur, 
e.g. for [// J = 1 the dipole oscillations are undamped 
for all N, see Figs. [Ha)-(b). The dipole dynamics is 
then studied as a function of the ratio U/J. In partic- 
ular. Fig. \3la) shows the damping rate F of the dipole 
oscillations as a function oi U/J for N = 11, 15 and 28 
[exemplifying cases a) , b) and c) above] . Here, F is calcu- 
lated using the expression for underdamped oscillations 
Xcomit) = e~^*[l - cos(f2t -I- (/)o)] + 2/0, with F, (j)o and yo 
fitting parameters. Three key observations are in order, 
i) The damping rate increases exponentially with U/J 
for intermediate interaction strengths 2 < t//J < 6, a 
result which is not captured by mean-field, and is signif- 



FIG. 2: Center of mass position as a function of time for the cases of Fig. [2a-c) and U / J = 1. The displacements 5 are 
indicated in the figure. The critical displacement 5c equals (5c = 6 and 2 in panels (a-b) and (c), respectively. 



icantly larger than what predicted using phase-slip tech- 
niques, vaHd for U <1 [isl. [lit. ^0 Eventually for large 
enough interactions [U / J ^ 6) the oscillations are over- 
damped for all N . Wc find that for the cases iV = 15 
and 28, this ovcrdamping corresponds to the formation 
of a Mott-state and a cake- structure as in Fig. [TJb) and 
(c), respectively. In particular, for = 15 the particle 
localization occurs for J7/ J « 4, a value remarkably close 
to the superfluid/Mott-insulator quantum phase transi- 
tion in an homogeneous lattice at commensurate filling 
and zero current. That is, the results for 5 < Sc nat- 
urally interpolate between the finite-current dynamical 
instability and the zero-current quantum phase transi- 
tion [l3l. iii) Despite the Mott-formation for large N, for 
a given U/ J the damping F is actually larger for smaller 
N, such that for iV = 11 the dynamics is frozen already 
for U/J < 4. In the following wc show that this has 
crucial consequences for the interpretation of the results 
of Ref. 0] in the most interesting regime of interactions 

U/J r^A. 

In the experiment of Ref. 01, the decay of dipole os- 
cillations was studied as a function of the optical lattice 
depth Vo for a fixed displacement (5 = 8, finding damping 
already for weak lattices Vo/Eft > 0.5, with Eji the recoil 
energy. The experimental data are shown as black dots in 
Fig-IHb) as a function of Vo in the range 2< Vo/Ep^ < 5, 
where the use of Eq. ^ is justified lH, [lj|j correspond- 
ing to the interesting regime of interactions 3 < J7/J < 8. 
For Vo/Er = 3 and Vq/Er > 3 the value of the damp- 
ing rate F has been extracted using formulas appropri- 
ate for underdamped and overdamped motion, respec- 
tively The most interesting experimental finding 
shown in Fig.[3jb) is the measurement of an abrupt tran- 
sition from a weakly damped regime to an overdamped 
regime for a lattice depth Vq/Ej^ ~ 3, where the damping 
rate F of the dipole oscillations increases by more than 
an order of magnitude. The physical mechanism behind 
this apparent transition has proven elusive. 

In Fig. [3jb) the experimental results are compared to 
our numerical results for iV = 80 and 45, green dia- 
monds and red squares, respectively. The value = 80 
has been chosen since it corresponds to the number of 



particles in the central ID lattice of the array in the 
experiment, which is the most largely populated with 
(rij) > 1 for all U/J, as in Fig. [Ijc). Conversely, the 
case iV = 45 exemplifies case (b), with (rij) < 1 for 
U/J > 4. The figure shows a very good agreement be- 
tween the numerical and the experimental results in the 
entire region 2 < Vq/Er < 5 (3 < U/J < 8). However, 
the case = 80 slightly underestimates the damping 
around Vo/Er ~ 4, while the agreement for = 45 
is almost perfect. For Vq/Er > 5 all numerical results 
fall inside the experimental errorbars, however, the case 
A^ = 45 shows a strong damping, while the case A^ = 80 
falls in the middle of the experimental errorbars. The ex- 
planation of the results above stems from the observation 
that in the experiment Sc varies between Sc ^ 18 and 15 
for 3 < Vo/Eji < 5, and thus S < Sc for all lattice depths. 
Wc can then use the results for the model systems of 
Fig. [31Ja) to explain the experimental findings. That is: 
i) the transition observed experimentally at Vq/Eh ~ 3 
is actually a crossover, where the ID systems with the 
lowest number of particles tend to localize first, in agree- 
ment with the discussion of Fig.[3l^a). ii) For Vq/Eji > 5, 
the dynamics of particles in the ID systems with (uj) < 1 
{N = 45 in the simulations) is completely frozen, and the 
overall mobility of the cloud is due to residual oscillations 
in lattices with higher onsite density. This latter observa- 
tion is in agreement with the results of Ref. , where it 
is shown that for Vo/Ef/ > 5 the damping rate observed 
in the experiment is well reproduced by the results for 
A^ = 80. We notice that numerical results for A^ = 80 
consistent with ours have been recently reported in [lo| . 
however the focus here is on a comprehensive explanation 
of the experiment 0] . 



The different behaviors of F for A^ = 45 and 80 and 
C//J > 4 can be modeled as follow. In the low-density 
case with A^ = 45 the tendency to localization is ex- 
plained by noting that interactions broaden the spatial 
width of the atom cloud, until the onsite density falls be- 
low one [see also Fig. [IJ a-b)]. In this case, the low-energy 
physics maps into that of an extended cloud of non- 
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FIG. 3: (a) Numerical results for the damping rate F of the dipole oscillations vs U j J for a fixed displacement 5 = 1 < 5c, with 
(5c = 6 (f2/J = 0.05623) and — 11, 15 and 28 [cases (a-c) in the text]; (b) Damping rate F for the experiment of Ref. 0] vs 
IJ I J and the lattice depth Vo/-E_f{. The experimental data, and the numerical results for A'^ = 45,80 are the black dots, the red 
squares and the green diamonds, respectively. 



interacting fermions, with single-band Hamiltonian [13| 

<ij> j 

with Cj and fermionic operators. For large enough 
displacements 5, the fermions largely occupy localized 
modes of the single-particle spectrum discussed above, 
and the COM remains frozen. The dynamics of interact- 
ing particles at large density, e.g. = 80 in Fig. [S^b), 
can be modeled starting from the case of largest interac- 
tions U/J ^ 1, where the density profile has a cake- like 
structure, Fig.[TJc). This situation is well described by an 
extended fermionization model [9., J^, , where Eq. ^ 
is replaced by an effective Hamiltonian with two coupled 
Fermi bands separated by an energy U [l9l | 

H2it) = -J [c\c, + 2d\d, + V2{4d, + d\c,)] 

<i,3> 

+ + sitrc]c, + + <5(t)]2 + u}d]d,), (2) 
j 

with the operators Cj , c] and dj , referring to the lower 
and higher energy bands of width 4J and 8J, respec- 
tively. Oscillations in this limit are due to the dynamics 
of the (delocalized) dj-fermions of Eq. ^ in the higher- 
energy band, while Cj-fermions are frozen in a (band) 
insulator. Observing these residual oscillations thus cor- 
responds to probing the superfluidity of bosons with two- 
particles per site in a homo gen eous lattice, in a local- 
density-approximation sense [20|. This picture, valid for 
U/J ^ 1 1^, [l^, can be extended to gain a qualitative 
insight in the dependence of the dipole oscillations on 
interactions for 4 < J7/J < 10. In fact, neglecting the 
parabohc potential, in this regime the model of Eq. ([2|) 
suggests that the spectrum is continuum, since the gap U 
between the two Fermi bands is smaller than their total 
width. It is thus plausible that Bloch-like oscillations of 



the particles are here suppressed, and transport restored. 
However, for U > 12J the energy spectrum develops a 
gap again around 4J, and thus transport in the lower- 
energy band is inhibited. Residual current is then due 
to delocalized particles in the higher-energy band, as ex- 
plained above. We notice that this picture is consistent 
with our numerical findings for U/J > 5 in Fig.[3Ub). 

In conclusion, we have explained the experiment in 
in the most interesting regime of intermediate interac- 
tions. The very good agreement between experimental 
and tDMRG results demonstrates the latter as a unique 
tool for quantitative comparisons with cold gases experi- 
ments in the strongly correlated regime in one dimension. 
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Clark are gratefully acknowledged. This work was sup- 
ported by OLAQUI, NAMEQUAM, FWF, MURI, EU- 
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code released within the PwP project (www.dmrg.it). 
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